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Abstract 

In this paper we study the system of two falling balls in continuous time. We modell the system 
by a suspension flow over a two dimensional, hyperbolic base map. By detailed analysis of the 
geometry of the system we identify special periodic points and show that the ratio of certain periods 
in continuous time is Diophantine for almost every value of the mass parameter in an interval. Using 
results of Melbourne m) and our previous achievements [T] we conclude that for these values of the 
parameter the flow mixes faster than any polynomial. Even though the calculations are presented 
for the specific physical system, the method is quite general and can be applied to other suspension 
flows, too. 


Introduction 

One of the main motivations for studying the statistical properties of hyperbolic dynamical systems is 
related to applications in physics. Yet, models with direct physics relevance are typically complicated, and 
thus results concerning such systems are quite limited. A remarkable exception is the class of dispersing 
billiards, see [6] for a detailed description of their theory. 

The system of falling balls investigated in the present paper cannot be regarded as a small pertur¬ 
bation of a dispersing billiard, especially as far as the dynamics in continuous time is concerned. This 
model introduced by Wojtkowski in m describes the motion of n point masses along a vertical half-line 
under the action of gravity, which collide elastically with each other and the floor. This system can 
be considered both in discrete and in continuous time. The first results were on hyperbolicity, i.e. non¬ 
vanishing of all relevant Lyapunov exponents. In m Wojtkowski proved hyperbolicity in case the masses 
of the particles are strictly decreasing up the line. Later Simanyi weakened this assumption, in [16] he 
showed hyperbolicity, when the masses of the particles decrease (but not necessarily strictly) up the 
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line and there are at least two different masses. Aiming at finer chaotic properties, in [12] Wojtkowski 
and Liverani developed a general method to show ergodicity for Hamiltonian systems and this could 
be used to show ergodicity of the system of two falling balls, in case the lower ball is heavier. Earlier 
Chernov showed a similar result in |3|. However, for three or more particles ergodicity is still an open 
question. In the ergodic regime of the two falling balls a detailed geometric description of the system and 
a quantitative analysis of the discrete time map lead to the verification of fine statistical properties in 
[Tj, in particular polynomial decay of correlations and the central limit theorem for Holder observables. 

In this paper we investigate the system of two falling balls in continuous time. Studying statisti¬ 
cal properties of hyperbolic flows is a technically involved task which is mostly related to the lack of 
hyperbolic behavior in the flow direction. In the past two decades there has been substantial activity 
in this field, here we summarize some results that are closely related to our work. An essential break¬ 
through was initiated by Chernov in [5|, where, using Markov approximations, he obtained a stretched 
exponential bound on time correlations for 3-dimensional Anosov flows that verify the so called ^uniform 
nonintegrability of foliations' condition (UNI for short). Dolgopyat developed Chernov’s result in two 
different directions. On the one hand, in |S] he showed that Anosov flows satisfying the UNI condition 
and a high degree of regularity are exponentially mixing. In his later work he studied the more general 
class of suspension flows over subshifts of finite type and in |7] he proved that such flows are typically 
rapid mixing, meaning that time correlations for sufficiently regular observables decay faster than any 
polynomial. Here the typicality condition is related to the presence of two periodic orbits such that 
the ratio of their periods satisfies a Diophantine condition. This condition plays an important role in 
our paper, thus we introduce the abbreviation DPO (Diophantine periodic orbits) for later reference. 
Dolgopyat’s result on rapid mixing was extended by Melbourne ([E]) to suspensions over Gibbs-Markov 
maps and also to suspensions over hyperbolic maps that can be modelled by a Young-tower (|IH|,[Ilj) 
with exponential tails. 

The DPO condition is a much weaker form of non-integrability in the flow direction than the UNI 
condition. Yet, in most applications, the verification of either UNI or DPO is based on the invariance 
of a canonical contact form or a perturbation thereof (see [TT| or [2] for example) Q Other than that, 
we are only aware of results that prove DPO - and hence rapid mixing - for a class of flows that is 
residual in an appropriate topology, and not for specific examples. Note that for Hamiltonian flows on 
cotangent bundles the canonical contact form associated to the symplectic form is preserved only if the 
Hamiltonian is a homogeneous function of the momenta (see cni, section 5.6). For the flow of two falling 
balls by the presence of a non-infinitesimal external field the Hamiltonian cannot be regarded as a small 
perturbation of a homogeneous function, hence we seek for alternative methods. 

In the present paper we prove DPO for the system of two falling balls, for almost every value of the 
mass ratio within a large interval (in the ergodic regime). To conclude that the system mixes rapidly in 
continuous time, we rely on |13j . This requires some additional work, as the periodic points originally 
constructed are a macroscopic distance apart, while for m it is essential that they are present on the 
base of the same Young tower. To obtain periodic points that can be realized on the base of the same 
Young tower, we apply a shadowing type argument, the details of which require most of the technical 
work in this paper. Similar ideas have already appeared in the literature, see in particular the notion of 
“periodic points with good asymptotics” in [S|. Nonetheless, in [S| good asymptotics is used to conclude 
stability of mixing and rapid mixing in the C” topology, while here we implement a shadowing type 
argument for the specific system of falling balls. 

As a consequence of our analysis we conclude that the system of two falling balls mixes rapidly in 


continuous time, for a set of mass ratios that has positive Lebesgue measure (cf. Theorem 2.2). It is 


worth pointing out that our analysis applies for almost every value of the mass ratio, as long as the 
mass of the lower ball is at least one and a half times larger than the mass of the upper ball. The only 
reason why we have to restrict to a smaller set of mass ratios is that in ^ the presence of a Young tower 
is established only for a smaller, yet open set of mass ratios. See our remarks after Formula (4.1) for 


^In that respect, dispersing billiard flows may be regarded as singular geodesic flows. 
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further discussion. 

The rest of this paper is organized as follows. In section we summarize the necessary prerequisites 
concerning the system of two falling falling balls, mostly from [T]. In section we discuss suspension 
flows and state our main results. Section contains the core argument of the paper, the construction of 
the periodic points satisfying DPO, along with the shadowing type argument, for mass ratios m G (|, |)- 
Finally, section discusses the extension to other values of the mass ratio. 


1 Setup and notations 

In this section we introduce the system and recall the necessary notations and results from our earlier 
paper [1]. The exposition is self contained, for further details about the dynamics and its properties we 
refer to our previous work. 

The system of two falling balls, introduced by Wojtkowski in describes the motion of two point 
particles of masses mi and m 2 that move along the vertical half-line, subject to constant gravitational 
force, and collide elastically with each other and the floor. We consider the case when the lower ball is 
heavier (i.e. mi > m 2 ), which corresponds to ergodic and hyperbolic dynamics (as shown in |12j and 
m)- As the action of ball to ball collisions depends only on the ratio of the two masses we rescale 
these masses such that mi -I- m 2 = 1. We introduce our mass parameter m and from now on we use the 
notation that mi = m, m 2 = 1 — m, where m € (1/2, 1) since we are in the ergodic case. 


gravity 


V2 



m2 


mi 


Figure 1: The system of two falling balls 

We neglect air resistance and assume all collisions to be totally elastic, therefore the flow preserves 
the total energy of the system, which we set to be 1/2 for practical reasons. We discretize time by 
considering the outgoing Poincare section corresponding to moments when the lower ball hits the floor 
and the next collision will occur between the two balls (and not between the lower one and the floor). 
Based on the work of Wojtkowski we use the coordinates 

h = -mvi and z = V 2 — vi, 

to describe the system, where Vi is the signed velocity of the Fth ball. This means that h is the energy of 
the lower ball (since it is on the floor it only has kinetic energy) and z is the difference of the velocities. 
Note that these coordinates are invariant between collisions. The phase space is then 
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Ml := |(/i, z) e M^l 0 < /i < 1/2, 1/2 — k > ^(1 — m) \ > 

to(1 — ~ — 2h + m < o|, 


where the conditions arise as follows. 

1. The first condition says that the energy of the lower ball is positive, but not greater than the total 
energy of the system, which we set to 1 /2 previously. 

2. The second condition is the inequality that implies that the upper ball has positive height. 

3. Finally, the third condition is to ensure that the two balls will collide before the lower ball returns 
to the floor. 


Recall that our Poincare section M\ corresponds to situations in which 

1. the lower ball is on the floor and 

2. it will collide with the upper one before returning to the floor. 

Hence, starting from a configuration {h, z) G Mi, first the two balls collide, and then the lower one will 
hit the floor several times before getting back to A^i. Let us denote by R{h, z) the number of bumps of 
the lower ball on the floor before returning to starting from the configuration [h, z). Then for any 
n G N we introduce 

Rn := {{h,z) e Mi\R{h,z) = n}. 

It is shown in [T] that none of these sets are empty. They are of course disjoint, moreover, even the 
closures of any two of them are disjoint provided their indices differ by more than 1. This way, in fact, 
the domains of continuity for the dynamics are identified, which provide a partition of the phase space 
Ml. 



Figure 2: The phase space of the dynamics 

Indeed, let us denote the dynamics of the system by T : Mi —>■ Mi- Then using the classical 
Newtonian laws of mechanics one can calculate that 
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T{h, z) = (mF, —(2n + 2)V^ — z) for (/i, z) S R„ 


(1.1) 


where 


F = 1 — h/m + az^, and a = 1 — 3m + 2m^. (1-2) 

The latter quantities often occur in the formulas, that is why we introduced extra notation on them. It 


can be seen from (1.1) that T is continuous, moreover on each F„, but the crucial dependence of T 

dRn 


on n shows that T is discontinuous on U()Tg 


The discontinuities occur at the curves 


rn := dRnCi dRn+i- (1.3) 

These correspond to configurations starting from which the two balls collide, then the lower one hits 
the floor n times and finally they land on the ground at the same time, on top of each other. Hence it 
is not clear whether the lower one reached the floor before their collision, or if it was the other way 
around. The two possible cases correspond to two different limits, one where the initial point (h, z) G r^ 
is approached from inside Rn+i and the other where it is approached from inside i?„. 

As we pointed out T maps each F„ diffeomorphically onto its image. The Jacobian of the dynamics is 

( —1 2maz \ 

V2{n+1) _ {2n+2)V2az I • (1-4) 

mVF VF / 

An important consequence of this formula is that det{DT{h, z)) = 1 and hence the normalized Lebesgue 
measure on Ali is an absolutely continuous invariant probability measure of the dynamics, which, by 
ergodicity, is unique. It can also be derived that DT{h, z) is a hyperbolic matrix at every point and that 
it contains a rotation by 180 degrees. 

Figure demonstrates how T maps the set onto its image. For further details about the regularity 
properties of T [1], section 3 is referred. The important fact that we will use in this paper is that T is 
uniformly hyperbolic (proved in [T] subsection 3.3), in the sense that there exists a forward invariant 
unstable, and a backward invariant stable cone field (C" and C| respectively) and these cone fields are 
uniformly transversal to each other. Curves 7 such that the tangent line T^"/ lies in the unstable cone 
C“ for every cc € 7 , are referred to as unstable curves. Stable curves are defined in an analogous way. 
Again in [1] it is shown that stable curves are increasing, while unstable curves are decreasing in the 
{h,z) coordinates. 




Figure 3: The action of the dynamics on F„ 

Concerning terminology it is important to distinguish stable manifolds as special stable curves. The 
stable manifold of a point x S A^i is defined as the curve IF®(x) such that for y G IF®(x) we have 
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d(T'^x,T""y) 0 exponentially as n —>■ oo. Equivalently, is a stable manifold if T'^W^ is a smooth 

stable curve for any n > 0. By the general theory of hyperbolic systems with singularities, a unique 
stable manifold of positive length passes through almost every x G Aii (see eg. [5] and references 
therein). Unstable manifolds are special unstable curves with analogous properties. 

In our previous paper we proved the following statements (cf. [T], subsection 1.2 and section 5.) 

Theorem 1.1. There exists an open interval I C (1/2,1) such that for any mass ratio m G I the discrete 
time map T : Aii ^ Mi can be modelled by a Young-tower with exponential tails. 

Theorem 1.2. If the system has subexponential complexity for some m G (1/2,1), then the discrete 
time map T : Mi —t Mi can be modelled by a Young-tower with exponential tails for this m G (1/2,1). 

We also recall from our previous work that there is an involution for the map T, i.e. there is a smooth 
map I : Mi —t Mi such that 


T-^ =IoToI, (1.5) 

on every smoothness component of T~^. This corresponds to the natural time reflection in the continuous 
time system and its action in the coordinates {h, z) is given by 


I{h, z) = (to(1 — h/m + az'^), z) = (mF, z). 


( 1 . 6 ) 


It can be deduced that I maps each Rn to T{Rn), and it follows from (1.51 that I maps stable curves 
into unstable curves and vice versa, moreover it maps stable manifolds into unstable manifolds and vice 
versa. Also, the curves /(r„) are the singularities of the inverse dynamics. Finally, let us recall from 
[T] the notations for the corner points of the sets and T(i?„). These are the intersections of the 
singularities r„, or the inverse singularities /(r„) with the boundary of the phase space. All the curves 
r„ are stable curves, in particular they are increasing and hence it makes sense to talk about their left 
and right endpoints. They are given, respectively, by the formulas 


Bxn{m) = {Bxhn{m), BxZn{m)) 

/m(—2m(n + 2) + 2n + 3)^ 


7 Z 2 


V 2(1 — m)n{n + 2) + 2 ’ y^l — (m — l)n{n + 2) 

A„(to) = {Xhn{m),Xzn{m)) 

to( 3 + 2n — 2m(n + 1))^ n+1 


(1.7) 


( 


2(n + 2)2 - 2m(n + l)(n + 3) ’ + 2)2 - m{n + l)(n + 3) ^ ’ 


The curves /(r„) are decreasing and so their left and right (or in this case rather top and bottom) 
endpoints are well defined, too. They are given, respectively, by the formulas 


Ixn{m) = {Ixhn{m),IxZn{m)) 
/ 771 


71+1 


V2(7i + 2)2 - 2m(7i+l)(n + 3)’ ^(n + 2)2 - 711(71 + l)(n + 3)^’ 

Ibxn{rn) = {Ibxhn{m),IbxZn{m)) 

771 71 + 2 


( 1 . 8 ) 


( 


2(1 — 771)71(71 + 2 ) + 2 — (777 — 1)71(71 + 2 ) ^ 
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Figure 4: The corners of the sets and T{Rn) 


2 Statement of results 

To state our results we model the continuous time dynamics of the system of two falling balls by a 
suspension flow over the uniformly hyperbolic base map T, discussed in the previous section. Using the 
classical laws of Newtonian mechanics one can calculate how much time is needed for the flow to return 
to the Poincare section Mi starting from the point (h,z) € Mi - It is given as 


T :Mi^ K+, 


T{h,z) 


(2n + l)v^ + 



2(to 


l)z for (h, z) G Rn- 


( 2 . 1 ) 


It can be shown that r is piecewise with the same discontinuities as the discrete time dynamics T 
(see [1], subsection 3.7 for details). The flow is then isomorphic to the following suspension. 

The phase space for the suspension flow is defined to be 


M = {{x,t)\x € Mi,0 < t < t{x)}, (2.2) 

and, after setting the equivalence relation {x,t{x)) ^ (T(a:),0), the continuous time action of the dy¬ 
namics is given by 


: Ai —>■ AI, ^*{x,s) = {x,s+ t)/. (2.3) 

Finally, the normalized Lebesgue measure on M is an ergodic invariant measure for this suspension flow. 

We define a class of observables in the following way. For a function u : A4 —>■ K we set its 77 -norm to 
be 

IklU := Ikiloo -Hsup|u(a;,'u) - v{y,u)\/d{x,y)’^. 

x^y 

Then for m G N and ?? > 0, let C^’^{M) the collection of u : TM —)■ M such that ||u||to,j 7 := + 

+ ■ ■ ■ + < CO; where dt denotes differentiation in the flow direction. 
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Definition 2.1. The suspension flow is rapid mixing if for any n> 1 there exist m > 1, C > 0 such 
that 

< C\\v\\m,r,\\w\\m^r,t~'^, 

M MM 


' • w o dLeb — v dLeb ■ 


I dLeb 


for every v,w € and t>0. 

Here we state our results. 


Theorem 2.2. The continuous time dynamics of the system of two falling balls is rapid mixing for 
almost every m £ I, where I C (1/2,1) is the interval from Theorem\l.l\ 


Theorem 2.3. Assuming subexponential complexity for the base map T (so that Theorem 1.2 can be 


applied), the continuous time dynamics of the system of two falling balls is rapid mixing for almost every 
me [2/3,1). 


3 Periodic orbits 

3.1 Outline of our strategy 

As already mentioned in the introduction, the core of our argument concerns the existence of periodic 
orbits with sufficiently diverse periods in the following sense. If a point a; S Afi is periodic for T 
with discrete period k (ie. T^x = x) then it is also periodic for the suspension flow with flow period 
Tk{x) = t{x) + t{Tx) + ... + t{T^~^x). Given two periodic points x and y with the same discrete period 
k, their (flow) period ratio is defined as Furthermore, recall (eg. from [3]) that a number w G K is 

Diophantine if it is badly approximable by rationals, ie. if there exist K > 0 and /3 > 1 such that for any 
p,q G Z, q ^ 0 we have \ujq — p\ > K\q\~^. The set of Diophantine numbers is of full Lebesgue measure. 
We will say that two periodic points have a Diophantine period ratio if is Diophantine. 

In his work [7] , Dolgopyat showed that mixing suspension flows over subshifts of finite type are rapid 
mixing if there exist two periodic orbits for the flow with Diophantine period ratio. Later in |13j Mel¬ 
bourne extended the approach to a large class of nonuniformly hyperbolic flows using operator renewal 
theory (see also m)- These flows are the continuous time analogue of the nonuniformly hyperbolic maps 
studied in [18j . He showed that if the system is not rapid mixing, then there must be some resonance in 
the roof function. In particular if there exist four periodic points with periods satisfying a Diophantine 
type relation, then the suspension flow is rapid mixing (see [T3] Theorem 2.6). The proof of this result 
contains ’’double inducing”, first reducing from the flow to a nonuniformly hyperbolic diffeomorphism 
and then from this to a uniformly hyperbolic one. However one could do this in one go, immediately 
reducing to a uniformly hyperbolic base map at the cost of having a larger roof function. The advantage 
is that by having only one inducing, there is only one complex variable to appear in the proofs (in)- 
Therefore, if one can find two periodic orbits on the base of the Young tower with Diophantine period 
ratio, then one can conclude rapid mixing for the flow (cf. |13j . Remark 2.7). 

We recall the formulas for the induced dynamics T : Af i —)■ A^i and for the roof function r : Af i —)■ M. 
They are given by 


T{h,z) 


T{h, z) = {mF,-{2n + 2)v^ - z) 


= {2n+l)V^ + 



2(m 


l)z 


for {h, z) G Rn 
for {h, z) G Rn 


(3.1) 


where F = 1 — h/m + az"^ and a = 1 — 3m -|- 2m^. We emphasize their dependence on the parameter m. 

When searching for periodic points the easiest attempt is to look for fixed points. Simply by solving 
the equation T{h,z) = {h,z) we get the candidates 







(3.2) 


F„(m) = 


-(n + 1) 


2(1 - a{n + 1)2) ’ ^l-a(„ + l)2 


Even though these are solutions for the fixed point equation, not all of them are physical solutions, i.e. 
some of the F„(to)’s may lie outside of the phase space and hence do not correspond to a valid physical 
configuration. However one can check that Eb(m) and Fi{m) are always physical solutions. Indeed both 
relations Eb)™) e i?o and Fi{m) G i?i are satisfied for every m G (1/2,1). We calculate the periods of 
these points in continuous time 


r(Eo(m),m) = 


2m 


yjim — 2m2 
and consider the ratio of the two periods 

T{Fi{m),m) 


T{Fi{m),m) = 


2m 


•^—3/4 + 3m — 2m2 


\/3m — 2m2 


(3.3) 


(3.4) 


T{Fo{m),m) ^-3/4 + 3m-2m2' 

This function is on [1/2,1], strictly decreasing on [1/2, 3/4] and strictly increasing on [3/4,1]. As 
a consequence the ratio (3.4| is Diophantine for almost every m G (1/2,1). Yet, to conclude that the 


continuous time system mixes rapidly we can not just directly apply the results of Melbourne, because 
the points Fq and Ej have a macroscopic distance in the phase space and so they might not be represented 
on the base of the same tower as required in [13) . 

Our strategy will be as follows. With the help of the natural partition of the phase space we switch 
to a symbolic space. Using the geometric properties of the map and the partition elements, for certain 
values of the mass parameter we guarantee the existence of a subsystem, which is a two-sided full shift 
on the two symbols 0 and 1. We construct a sequence of periodic points {P„(to)} that accumulates on 
Foim), but in a way that the trajectory of P„(to) spends more and more time in the vicinity of Ei(m) as 
n increases. We also show that all points FQ{m) and {P„(m)} have sufficiently long stable and unstable 
manifolds, hence a tower can be built such that for large enough n both Foim) and Pn{m) are on the 
base of it. In terms of m. for rapid mixing it is then enough to show that the ratio of the continuous 
periods of P„(to) and Po(m) is Diophantine. 


3.2 The full shift as a symbolic subsystem 

For a point P G Ali we adjust the two-sided infinite sequence x = {xi} as its symbolic representation, 
where Xi = n iS T^{P) G P„. We denote the natural projection from the symbolic space to A4i by tt, 
set the past and future separation time for two sequences x and y as 


s-{x,y) = minjlfc] : A; < 0,Xfc ^ yk}, 
s+{x,y) = minj/c ■■ k>0,Xk yk}, 

and also the separation time as 

s{x,y) = minlJA:] : Xk 7 ^ yk] = min{s_(x, j/), s+(a:, y)}. 

We define the symbolic distance as 

dsymix,y) = e^^^’y), 

where 6 G (0,1) is chosen in such a way that the distance on Ali and in the symbolic space are related 
as 

d{TT(x),Tr{y)) < C ■ dsymix,y), (3.5) 
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Figure 5: The geometry of the sets Ri n T{Rj) for i,j G {0,1}. 


for some constant C > 0. By uniform hyperbolicity of the discrete time system (proved in [T]) there 
exists such a 9 G (0,1). 

We now show the existence of the two-sided infinite, full shift subsystem. 

Lemma 3.1. For every m G [2/3,3/4] the projection of any sequence x G {0,1}^ is realised as a physical 
configuration, i.e. as a point in M.\. Any such point has local stable and unstable manifolds that fully 
cross the phase space. 


Proof We will call a region quadrangular if it is diffeomorphic to a square and we will refer to such a 
region as a (curvilinear) rectangle if it is bounded by two stable (increasing) and two unstable (decreasing) 
curves. These will be referred to as the stable and the unstable sides of the rectangle, respectively. 

Consider the geometry of the sets Ri C T{Rj) for i,j G {0,1}. While Rq C T{Rq) is quadrangular 
for every m G [1/2,1) this does not necessarily hold for the other three sets. They can be triangular, 
quadrangular or pentagonal (in the previous topological sense) depending on the value of the mass 
parameter m. The different possible cases are shown on Figure]^ We will work only with the case 5(c)[ 
when 2/3 < 771 and hence all four regions are quadrangular. The sets of points that correspond to the 
symbolic representations {x : xq = 0} and {a; : ccq = 1} are of course Rq and Ri. We go one step 
further and identify the sets corresponding to symbolic representations {a; : a:o,a;i G {0,1}}. Formally 
they are Ri C T~^(Rj) for i,j G (0,1}, i.e. the preimages of the four quadrangular regions we have just 
discussed, but what is important is again the geometry of these sets, see Figure]^ They form (curvilinear) 
subrectangles within Rq or i?i, fully crossing them in the stable direction. 

Actually, in what follows we prove a somewhat stronger claim. 

Claim The stable sides of the rectangles Ri C T~^[Rj) for i,j G (0,1} cross both unstable sides of 
both of the rectangles T{Rq) and T(i?i). 

In fact this holds automatically in case of the regions i?i C T~^{Rq) and Ri C r“^(i?i), because 
their stable sides are the preimages of segments of certain singularity curves - xq, ri and the segment 
of dM .1 that forms the top edge of T{Ri) - which connect the two unstable sides of T{Ri). Under the 
action of the inverse dynamics these unstable sides map onto the unstable sides of Ri , which are bits of 
the boundary of the phase space. Consequently, the above mentioned segments of the singularities map 
onto stable curves stretching from side to side in A4i and crossing the unstable sides of both T(i?o) and 
T(i?i). 

The situation for the regions i?onT ^(i?o) and i?onr ^(i?i) is slightly more complicated, essentially 
because Rq is topologically a triangle and so we can not talk about its stable and unstable sides. Actually 
it may happen that the top edge of Rq C T“^(iii) does not intersect the left side of T{Ri) as required. 
To identify the cases when this happens we calculate the left endpoint of this top edge and compare it 
to the top left corner of T{Ri). As follows from our discussion in sectionFormulas (1.7), (1.8) and 
(1.1), the upper left corner of Rq C T“^(i?i) is the point T“^(Ai(r? 7 )), its second coordinate is 


The upper left corner of T{Ri) is the point Ixi{m), it has second coordinate ^ 7 = 


8m ' 


8m ' 

It follows then 
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that for TO < 3/4 the point T~^{Xi{m)) falls outside the region T{Ri). Since Rq n T~^{Ri) lies above 
the region Rq n T“^(i?o) this implies that for 2/3 < to < 3/4 the stable sides of both regions cross the 
unstable sides of both T{Rq) and T(i?i). This completes the proof of the Claim. 

If we want to proceed one more digit and identify the set of points corresponding to sequences 
with the first three digits arbitrarily chosen from {0,1}, we should consider the preimages of the sets 
T{Ri) n Rj n T~^{Rk) for i,j,k G {0,1}. Note that due to our previous argument these sets form 
subrectangles in T(i?o) and T{Ri) fully crossing them in the stable direction, but narrower in the 
unstable direction (see again Figure [^. This, together with the previous geometric observations, implies 
that the preimage of any of these subrectangles is a subrectangle in one of the regions Ri C T~^{Rj) for 
i,j € {0,1} fully crossing it in the stable direction. It follows that the process can be iterated showing 
that for any to G [2/3,3/4] and any given one-sided infinite sequence Xo,Xi, ■ ■ ■ G {0,1} the projection 
of the set {y ■ yi = Xiii G N} is not empty. It is actually a curve, moreover a fairly long local stable 
manifold, that crosses the phase space. 

To complete the proof it is enough to take into account that the involution / : AIi —>■ AIi maps the 
set F T(i?j) to itself and stable manifolds into unstable ones. The action of / (more precisely 

the action of it lifted up) on the symbolic space is given by 

/(x)i = 

Hence for any sequence x G {0,1}^ the projection 7r(3;) is realised as a physical configuration. More than 
that it has long local stable and unstable manifolds 

W%tt{x)) = TT{{y :yi=XiyiG N}), W'"{tt{x)) = TT{{y : yi = x^Wi < 0}), 

crossing the phase space. □ 


II 





3.3 Convergence of period ratios 

With the help of this symbolic subsystem constructed in the previous subsection we can find the periodic 
points needed for our purposes. First of all consider the two special sequences, one consisting of all zeroes 
the other of all ones. Let us denote them by 0 and 1, respectively. Their projections are 


= Fo{m) 7r(l) = 


(3.6) 

where the Fi{mys are the fixed points identified in (3.2). We now define our sequence of periodic points 
for any m S [2/3, 3/4] as 


P„(m) := 7r(p") ,where p/ = 0 for — n < i < n — 1, 
p/ = 1 for n < z < n + — 1 

and p is periodic otherwise. 

Defined in this way Pn(rn) has discrete time period + 2n and since 

s(Fo(m),P„(m)) = n 


(3.7) 


it is exponentially close to Fo(to) by Formula (3.5). The local stable manifold of Pn{rn) intersects the 
local unstable manifold of Fo{m) and the same holds with stable and unstable roles exchanged (this is 
actually true for any two points with symbolic representations from {0,1}^). From this, together with 
the second half of Lemma |3.1[ it follows that for n large enough the points P„(m) and Fq (to) can be 
represented on the base of the same Young-tower. What remains is to show that in terms of period in 
continuous time P„(to) behaves more and more like Pi (to) as n increases. We make this precise in the 
next two propositions. 

Proposition 3.2. Let denote the k-th Birkhojf sum of t. Then for any to S [2/3, 3/4] we have 

Tn^+2niPnirn),m) T{Fi{m),m) 

(n'^ + 2n)T{Fo{m),m) T{FQ{m),m) 

as n ^ oo and the convergence is uniform in to. 


Proof. First we show that r is uniformly bounded on the set (i?o U i?i) n (r(i?o) U T{Ri)), cf. (3.8). 
Recall the formula for r given in (3.1). On the phase space we have the following trivial bounds 


h<l/2 z<0. 


The relation between the first coordinate of the involution and the quantity F (see (1.6)), together with 
the previous bound on h implies that F < 1/(2to). Finally on the set we are working on n = 0 or 1. 
Substituting the previous bounds into (0 gives that 


t(/i, z, to) < < 


4v^ 

75 ' 


(3.8) 


We want to relate the value of the roof function r along the trajectory of P„(to) to r(Pi( to),to). 
By construction the trajectory of P„(to) gets exponentially close to Pi (to) and stays close to this fixed 
point for an order of iterations. Hence to compare the two values we are going to perform a Hrst order 
expansion of t in the space variables (not in the parameter to). 

The relevant partial derivatives are 


dr 

dz 


1 


'/2/rnh 

= 2(1 — to ) 


(2n + 1) 
+ (2n -k 1) 


-1 


m\/2F 

2az 

TP’ 


(3.9) 
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where again F and a are as in (1.2 1 . To estimate these quantities on the set (i?oUi?i)n(T(i?o)UT(i?i)) we 
use again that n = 0 or 1. Furthermore, we need lower bounds on h and F (because they appear in some 
denominators) and also on z (because a < 0 for m G [1/2,1]). Now a lower bound on h can be the first 
coordinate of the leftmost point of our domain, which is Ixhi{m) = ■ Using again the relation 

in (1.6), together with the fact that our domain is mapped onto itself by I, we get that F > 2 ( 9 - 8 m) • 
Finally the minimum of z is given by the second coordinate of the lowest point of the domain. Since 
this is a point of ri, which is an increasing curve, it can be further estimated by the second coordinate 
of the left endpoint of ri, which is Bxzi (m) = ^ 4 ! 3 ^ • It is then clear that all estimates are continuous 
functions of m on [2/3, 3/4] (actually on the whole parameter domain [1/2,1]) and all denominators in 
(3.9) are separated from 0, hence there is a uniform bound on ||V/i_ 2 ;r||. Let us denote this bound by 




Now by the mean value theorem the continuous time period of Pn{m) can be expressed as 


r„2+2n{Pn{m),m) =Tn{Pn{m),m) +T„(T”"+"P„(m),m) + n^r(F’i(m), m) + 

n^+n-l ( 310 ) 


Here Uj is the unit vector parallel to the line connecting F’i(to) and T'^Pn{m), Dy denotes differentiating 
in the v^ direction and is the point on the segment specified by the mean value theorem. Apart from 


the bounds (3.8) and those given on the gradient of t we will also estimate d(Fi(m), T*P„(to)) using 
(3.5). For this note that 


s(Fi(m), T*P„(m)) = min{i — n,n'^ + n — i} for n < i < + n. 


so the sum of the distances in (3.10) can be estimated by segments of two convergent geometric series 
(one for the indices n < i < nrj2 + n and the other one for n^/2 + n < i < + n — 1), both with 

quotient 9. Hence we have the overall estimate 


Tn2+2niPn{m),m) _ r(Pi(m),m) ^ 

(n^ + 2n)r(F'o(w),m) T{Fo{m),m) ~ 

^ n^T(fi(w),TO) _ T{Fi{m),m) ^ 

~ (n^ + 2n)T(Fo(m),m) T{Fo{m),m) 

2n-\s/W^ G'oo2C(l - 6>"U^+^) 

(n^ + 2 u)r(Po(w), m) (n^ + 2 n)r(Po(”^)j w)(l — 0 ) ’ 

from which the statement of the proposition immediately follows. □ 

To complete our argument we have to show that not only the ratio of the two periods converge but 
also its derivative with respect to the parameter m, as expressed in the following Proposition. 

Proposition 3.3. For any m G [2/3,3/4] we have 


d j 

^ T„2+2n(-P„(m),m) / 

1 ^ 1 

( T(Pi(m),m)\ 

dm ' 

^ (n^ + 2n)T{Fo{m), m)) 

^ dm' 

/t(Po(w),to)/ 


as n ^ 00 and the convergence is uniform in m. 

The proof of Proposition |3.3| requires more work, hence we move it to subsection |3.4[ 
Now we are in a position to prove our first result. 
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Proof of Theorem \2.2\ We have the sequence of special periodic points converging to Fo{m) for 

every m G [2/3,3/4]. As stated in Theorem in [T] we showed the existence of an open interval in 
(1/2,3/4) such that for every m G / the discrete time map T : M.i M.\ can be modelled by a Young 
tower with exponential tails. The size of such a tower is determined by certain constants of the system, 
for example uniform bounds on the curvatures of unstable manifolds and of singularity curves, and also 
the minimum rates of expansion and contraction. What is important is that for m G [2/3,3/4] these 
constants can be chosen independently of to, since the dynamics is and the dependence on to is 
continuous. Then we have a tower with uniform size, hence for n large enough Pn{m) and FQ^m) are 
close enough to each other and so, as we have shown, they can be represented on the same Young-tower 
for any value of to. Applying the result of Melbourne m) it is then enough to prove that the ratio of 
the continuous time periods of Pn im) and Fo(to) (the latter considered under (n^ -|- 2n) iterations of T) 
is Diophantine. Propositions 3.2 and 3.3 imply that this ratio as a function of to converges in the 
topology to the function in (3.4). The limit is a C^, strictly decreasing function for to G [2/3, 3/4], hence 
by further increasing the value of n (if necessary) the ratio +2n{Pn{rn),rn)/{{n^ + 2n)T(Fo(TO), to)) 
will also be strictly decreasing and therefore indeed Diophantine for almost every to in this interval. □ 


3.4 Proof of Proposition |3.3| 

Differentiating the fractions appearing in the statement of the Proposition we get 


d f Tn-^+2n{Pn{m),rn) \ ^ 

dm \ (n^ + 2n)T(To(TO), to) J 

T(fo(TO),TO)g^T„2+2»(Pn(w),TO) - r„2+2n (-Pn (w) , to) g^T(To (to) , to) 
(rP + 2n)r2(Fo(TO), to) 


d / T{Fi{m),m) \ _ 
dm \T{Fo{m),m) J 

T{Fo{m),m)-£^T{Fi{m),m) - T{Fi{m),m)-£^T{Fo{m),m) 
r 2 (Fo(TO),TO) 


Hence, by Proposition 3.2 it is enough to prove that 


Proposition 3.4. 

2 I o -l—Tn^+2n{Pnim),m) -)■ — T (i^i (to) , to) , 
+ Zn dm dm 

as n —>■ oo. 


To show this we first derive a useful formula for the derivative. Before that let us introduce some 
notations for brevity. For any k G {0,..., -|- 2n — 1} let 

P^{m) :=r'=P„(TO), 

so that for any integer j we have Pn{m) = Pl^{m), where I = j {mod rP + 2n). Also we will refer to 
the two dimensional vector obtained by differentiating the coordinates of P^{m) with respect to to as 
^P/(to), and ^^P/(to),1 ^ will denote the three dimensional vector that has first two coordinates 
identical to ;^P/(w) and 1 as the third coordinate. Using this notation we can write 
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(3.11) 


-^Tn2+2n{Pn{'m),m) = ^ —r(P^(m), to) 


n^+2n-l 


n~ -\-zn — 1 j , 

^•—n \ 




where (.,.) denotes the usual scalar product in K^. While we have an explicit formula for r and hence 
also for '^h,z,mT: calculating the two dimensional vector ^P^(m) is slightly more difficult. For this we 


use that 


T” = P^{m), 


and we perform implicit differentiation with respect to m keeping in mind that the dynamics T also 
depends on m. Let us denote by §i^{Pn^^) the two dimensional vector obtained by first differentiating 
T with respect to m and then evaluating the result at the point As we perturb the parameter we 

have to take into account that not just every point of the trajectory varies with to, but the dynamics 
also chenges. These two effects jointly appear in the calculations leading to the formula 


-^Pnim) = 
dm ^ 


n^+2n—1 \ n'^+2n—1 n"^+2n—1 

= |/- n DT{P^+\m))\ Y. n DT{P^+^{m)) 

j=i+l 


i=0 


z =0 


dT 

dm 




= [I - DT 


m +2n 


(Puir 


n'^+2n—1 


Y DT” 


^+2n—(i+1) ^pk+i+1 


2=0 


n'^+2n—1 


8T 

{m))-—{P^+\m)) 


dT 


= Y +2-)(p„^(TO))-pr+i(p„^(TO)))-i.—( p„"+“(to)), (3.12) 


i=0 


where in the last line we used the inverse differentiation rule and again the fact that the point P^{m) is 
periodic with period + 2n. In the next step we show that this sum is bounded. 


Lemma 3.5. For each k G {0,..., + 2n — 1} the quantity g^P^ (ui) 


is uniformly bounded. 


Proof. To make the notations simpler in this proof we will suppress the dependence of the objects on 


the parameter to. Let us denote the i-th term in the sum (3.12) by Vi(n, k), so 


= (£, 7 ^d+l)-(«"+ 2 n)(pfe) _ dT^+\P^)) ■ v,{n,k). 


(3.13) 


Observe that because P^ is a periodic point the tangent spaces of its stable and unstable manifolds can 
be calculated as the stable and unstable eigendirections of the tangent map DT” ^^"(P,J), respectively. 
We denote the normalized eigenvectors of this matrix by s(n, k) and u{n, k) where the letters u an s 
refer to stable and unstable, respectively. Then by the invariance of these directions we have 

DT{P^)s{n, k) = /r(n, k)s{n, k + 1), DT{P^)u(n, k) = \{n, k)u{n, k + 1), 

defining the quantities X{n,k) < —1 < fj.(n,k) < 0, (both are negative as the tangent map contains a 
rotation by 180 degrees). Here and throughout the subsection the index k should be understood modulo 

n^+2n—1 n^-\-2n—l 

+ 2n. Note that ])[ fi{n, k) and H A(n, k) are the stable and unstable eigenvalues of the 
fe=o fc=o 
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matrix Z3T" and hence their product is 1, since I?T" has determinant 1. We consider 

the decompositions 


dT 

Vi(n, k) = ais(n, k) + huin, k), -—= Cis(n, k + i + 1) + d^uin, /c + i + 1), 

om 


and substitute them into (3.13) to get 


Cis{n, fc + J + 1) + diu{n, k + i + 1) = 

n^+2n—(i+1) ^ n^+2n—(z+l) ^ 

= “• n + ' + +'■■ n ipiyiTT) 


u(n, k + i + 1) — 


— Oi ^(n, k + j)s(n, fc + i + 1) — 6^ A(n, k + j)u{n, k + i + 1) = 


3=0 


3=0 


I /n +2n \ 

= //(n, k j)( PJ A(n, A:) — 1 j s{n, ki1) + 

j^o ^ k^l ^ 

i yn^+2n \ 

+ hY[X{n,k+j)( [] /i(n, k) — Ij u{n, fc + i + 1 ). 

w—n \ i —1 / 


This gives the relations between the coefficients 


Gi = 


Ci 


1 


= 


i n2+2n 

lld{n,k+j) n A(n,fc)-1 

J=0 k=l 

di 1 


(3.14) 


i n'^+2n 

U.Kn,k + j) n Kn,k)-l 

3=0 k=l 


First we show that both Ci and di are uniformly bounded. To see this, by the uniform transversality of 
stable and unstable cones it is enough to check that || f^(Pn'''*)|| is itself bounded. We have the general 
formula for the derivative of (3.11 

dT 


^{h,z,m) = (l + z^{l - 6m + 6m"), + z"(4to - 3))) , 




and as we have already shown right after (3.9), on the domain {Rq U Pi) n (T(Po) U T(Pi)) and for the 
parameter interval m G [2/3, 3/4] the quantity F is bounded away from 0, while h and is bounded, 
hence the whole norm ||f^(Pn^*)|| is bounded, too. 

Finally, since the dynamics is uniformly hyperbolic for all fixed values of m (proved in [1]), there are 
numbers fi and A such that A(n, fc)<A<—l</r< /r(n, fc) < 0 for all k. Hence both sets of numbers 
{a„ 2 _|_ 2 „_i_i}(Tg (i.e. the o^’s in reversed order) and can be estimated by the initial 

segment of a geometric series. Indeed this is immediate for the bi’s, and for the a^’s note that when 

i 

the index i is large and hence the denominator ]/[ fj.(n, k + j) is very small, it is still compensated by 

3=0 

n^+2n 

]/[ X{n,k). Therefore these sequences are summable and so elementary inequalities imply that the 
statement of the lemma is true. □ 
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Now we are in the position to prove Proposition |3. 4 


3.4 


Proof of Proposition\3.4\ C onsi der the formula (3.111 we gave for the derivative on the left hand side of 


By Lemma 3.5 we know that each vector {-£-Pf{m), 1) is bounded and we have already 


Proposition 

shown after (3.9) that the gradient of r is also bounded, however this was the gradient only in the spatial 


variables. In (3.11) the gradient contains the derivative with respect to m, too. By (3.1) 


dm 


2n + 1 




h 

2m3‘ 


-2z, 


and basically the same argument given after (3.9) shows that this quantity and hence the whole gradient 


in each of the terms of (3.11) is bounded. Therefore each scalar product in (3.11) is bounded as well. In 
(3.11) one can replace P„(to) by Fi{m) to obtain a similar formula for . After this 


dTn2^2n{Pnim),m) dT{Fi{m),m) 


T 277. 


dm 


n^+2n—1 


2n 


z=0 


dm 




dT{Fi{m),m) 

dm 


+ 2n 


E 

i— 3 n 


i\'h,z,TnT)iPnim),m) - (yh,z,niT){Fi{m),m), ( —P*(to),1 


+ 


n —n 


+ 2n 


i— 3 n 


(v„ 

,2,m 'r)iFi{m),m), 


^Pnim) - ^Fi{m),0 
am am 


0(1). (3.15) 


Note that we compressed an order of n number of summands into the error term and so we are 
left to deal only with the middle part of the original sum. The first sum is easy to handle. Basically 
it is enough to check that '^h,z,mT is in the variables {h,z), which turns out to be the case after 
doing similar computations as before. After this an argument similar to the one used in the proof of 
Proposition 3.2 shows that the first sum in (3.15) tends to 0 as n —>■ oo. Actually it is exponentially 
small in n, because for 3n < i < n? — n the sum of the distances between the points Pf{m) and Fi(m) 
can be estimated by a segment of a geometric series with the largest term being exponentially small, 
since s{Pf^{m), Fi{m)) > 2n for such indices. 

What remains is to show that the quantity 


1 


n —n 

2n E/ 

i— 3 n 


{'^h,z,mT){Fi{m),m), ( J-Pf{m) - ^Fi(m),0 
' dm dm 


(3.16) 


tends to 0 as 77 —)■ oo. We know that the gradient iyh,z,mT){Fi{m),m) is bounded so we estimate the 
norm of the other vector, which is a difference of two derivatives. We use the explicit formula (3.12) for 
these differentials to get 


^Puim) - = 

dm dm 


n +2n—1 


dT 


= ^^^HPnim)) - DT^+^{Pf{m)))-^ ■ ^(P^+*(777))- 


j=0 


- (£)r(^+i)-(”'+2")(Fi(777)) - Dy+\Fi{m)))-^ ■ E(Fi(r77)). (3.17) 
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To estimate the middle part of this sum we use ideas from the proof of Lemma 3.5 (especially formula 
(3.14)), which remains valid and actually is simpler, when P^(m) is replaced by Fi{m). From there 
it follows that for n<j<n'^+n—1 the terms that form the j-th difference in (3.17) are already 
exponentially small in n so we can omit them. Therefore it is enough to estimate the differences in (3.17) 
for ‘Hn < i < V? — n, 0 < j < n — \ and + n < j < n? + 2n — 1. Using the triangular inequality we 
estimate the norm of the j-th term by 






dT 

dm 


{P^^\m)) 


(3.18) 


Let us denote the vector in the first term by Note that ^ is on the domain we work on, so 
after setting 

dT dT 




we have 




< C ■ d{Pt\m),Fi{m)) < 


0 ) 


(3.19) 


We denote the coordinates of Vj^i and w,- 




in the basis {s(Pi(m)), u(Fi(m))} by {aj^i,bj^i) and {cj^i,dj^i), 


respectively. Then a calculation analogous to the one in the proof of Lemma 3.5 (c.f. (3.14)) gives 


_ Uj-y ■ A(fi(TO))-^+^ 
~ A(Fi(m))""+ 2 " - 1 

d - = ■A^(Pi(w))J+^ 

/r(Fi(m))”"+2’"-1' 


(3.20) 


Due to the uniform transversality of stable and unstable cones, estimates on the norm of Wj^i are - up to 
a constant - the same as estimates on the coordinates Cj^i, dj^i. By the same reason (3.19) implies that 

\aj,i\ < C ■ < c ■ 0s(Pi+'im),Fi(m))^ 


in particular they are bounded and the estimates form segments of geometric series in j. Hence the sum 
of Cj/s for 0 < j < n — 1 and the sum of dj/s for + n < j + 2 n — 1 are both exponentially 

small in n. To see that this holds also for the sum of Cjy’s for + n < j < + 2n — 1 and the sum of 

dj/s for 0 < j < n — 1 use in addition that 3n < i < n? — n and hence s{PP^'^{m), Fi{m)) > n. These 
observations altogether leads to the fact that for all i fixed between 3n and n? — n 


n—1 n^+2n—1 

I 

j—O j—n^-\-n 


wj4<c-0^. 


Therefore it remains only to estimate the second term in (3.181. For this we redefine the coefficients 

^j,i ? • • • 5 ^j,i 


— (py*(m)) = aj^is{P^{m)) + hj^,u{P^{mj) 

dT 

vj— (-Py*("^)) = Cj^is{Fi{m)) + dj^iu{Fi{m)). 


(3.21) 
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In this notation, after a calculation analogous to the repeatedly referred one from Lemma |3.5| the second 
term of (13.181) reads as 


1,'i+i+l 


j n^-\-'2Ti 

n A(p„^(m))-i 

fc=0 k—l 


s{Pn{m)) + 




j n^+2n 

nA(p^'=(m)) n Kp^{m))-i 


u{Pnim))- 


k—0 k—l 

_ ■ \{Fi{m)y+^ 

A(i^i(m))""+2" - 1 




(3.22) 


To finish the proof we use the regularity of the stable and unstable directions stated in the following 
lemma. 

Lemma 3.6. Let x and y be two points with both stable and unstable directions well-defined. Further 
assume that for all i for which —s-{x, y) < i < s+{x, y), the points T'‘(x) and T^(y) are either in Rq or 
in Ri- Then there exist constants Cs,Cu > 0 and 7s,7u G (0,1) such that 

||s(a:) - s(j/)|| < ||u(x) - u{y)\\ < 


This Lemma expresses the dynamical Holder continuity of the stable/unstable directions, but we 
were unable to find a good reference, so instead we give a proof of the Lemma in the Appendix. As a 
consequence, using also that ||ldT|| is bounded on our domain, this property holds for the functions /i(.) 
and A(.) as well. Finally note that all coefficients aj^i,... ,dj^i can be expressed using scalar products 
of the vectors u(.), s(.), their orthocomplements and We have already shown that the latter is 
uniformly bounded and this, together with the uniform transversality of stable and unstable cones, leads 
to the fact that both differences — cj^i and — dj i are exponentially small. It is then a 

straightforward calculation, using several triangular inequalities, that (3.22) is exponentially small and 
hence so is the second term in (3.18). This completes the proof of Propositions 3.4 and thus 3.3 □ 


4 Extension of results 

As the reader might have already noticed the calculations in Section do not depend on the system in 
hand crucially. The system specific parts of the argument were 

1 . to find a full shift subsystem based on the geometrical properties, 

2 . to give bounds on certain quantites like r and its derivatives, or the expansion and contraction 
rates of T, 

3. to show that the ratio of the continuous periods of the original fixed points for T is a function 
of the parameter m with nonzero derivative. 

Appart from these the whole method works in full generality. Taking advantage of this observation, in 
this section we extend our results for a larger set of parameters. 

As for 2., without going into details or doing the actual calculations, we claim that whenever we work on 
a domain of the form F T{Rj))-: £^11 the mentioned quantities will be uniformly bounded 

for the corresponding interval of the parameter. This can be checked in an analogous way as we did in 
Section [3] for k = 1. 

In Lemma [3. 1 1 we discussed the geometry of the sets Ri n T{Rj) for G {0,1} and showed that their 
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preimages intersect the unstable sides of the sets T{Ro) and T{Ri) if m G [2/3, 3/4]. From this we 
could conclude the existence of a subsystem that is a full shift on the two symbols 0 and 1. Now we 
can essentially repeat our argument for the stripes Rk-i,Rk with larger k giving a similar full shift 
subsystem for different parameter intervals. This is even simpler than in Lemma |3.1[ because for fc > 1 
all stripes Rk are topologically squares, in contrast to Rq that is topologically a triangle, which required 
special care. 

Lemma 4.1. Given any k > 1, for every m G ^ 2 fc^+ 2 fc^ ] projection of any sequence x G 

{k — 1 , is realised as a physical configuration, i.e. as a point in A4i. Any such point has local stable 
and unstable manifolds that fully cross the phase space. 


Proof. The k = 1 case was proved in Lemma |3T| To apply the same argument to the case k >2, one has 


to check only that the sets RiC\T{Rj) fori,/ G {k—l,k} are all quadrangular. The required intersections 
and hence the whole iteration scheme is then automatic. Also, checking that all the previous four sets are 
quadrangular is easy. The geometry of the system, in particular that the singularities are all increasing, 
while the inverse singularities are all decreasing curves, implies that it is enough to show that /(r^.) 
intersects both rk and rk- 2 - A nonempty intersection rk Gl{rk) can be guaranteed, using the mentioned 
monotonicity of the curves, by showing that the left endpoint of has smaller h coordinate than the 
right endpoint of I{rk). Based on (1.7) and (1.8) this condition is equivalent to 


Bxhk(rn) < Ibxhk{m) ^ 




m{—2m{k + 2 ) + 2 fc + 3)^ 
2(1 — m)k{k + 2 ) + 2 

1 + fc 

2 + k 


< 


m 




2 (l-TO)fc(fc + 2) + 2 


< m. 




For the intersection of rk -2 and I{rk) it suffices to check that the left endpoint of rk -2 has smaller z 
coordinate than the left endpoint of I{rk). Using again (1.7) and (1.8) this leads to the formula 


- 


Bxz{m, k — 2) < Ixz{m, k) 
k ^ _ k + l 


(m — l)k{k — 2) •\/(fc + 2)2 — m(fc + 1)(A: + 3) 

2 /c 2 +2k-l 






2 fc 2 + 2k 


> m. 


This completes the proof of the Lemma. 


□ 


As a final step towards the extension of results we need to present fixed points depending on to, such 
that their continuous period ratio is a function of to with nonzero derivative. Luckily the fixed points 

identified in (3.2) will do the job. Indeed, consider the fc’th parameter interval from Lemma 4.1 i.e. let to 
1 ^ 


[i+l’ ^ 2 fc^+ 2 fc ^ ] ■ Then the fixed points Fk-i{m) and Fk{m) are realised as physical configurations, 
they correspond to the projections 7 r( fc — 1 ) and 7 r(fc). The ratio of their periods is 


T{Fk{m),m) [k + l)Vl — k'^a 


T(Ffc_i(TO),TO) fci/l - (fc + l)2a’ 


(4.1) 


where a is as in ( |1.2[ ). Recall that in Theorem 2.2 we stated rapid mixing for almost every value of 
the parameter within the interval I C (1/2,1) from Theorem 0 the existence of which was shown in 
our earlier paper [1]. We do not have a quantitative description of I, we only know that it is an open 
interval containing 0.74 (though computer assisted techniques could give bounds on the endpoints of 
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I). In any case, the proof of [T] that the discrete time system can be modelled by a Young tower with 
exponential tails works only for parameter values within I. On the other hand, as stated in Theorem |1.2[ 
if we knew that the singularity set has subexponential complexity, then we would have a Young tower 
representation. Hence to extend our results we add the extra assumption of subexponential complexity (a 
usual assumption in the literature ensuring that expansion prevails cutting), which in particular implies 
the growth lemma and therefore the required tower representation. 


Proof of Theorem \2.S\ In view of our previous analysis the proof of Theorem |2.2| can be repeated for 
any of the parameter intervals [ 5 ^, ^ 2 fc^+ 2 fc^ ] on.ce we show that ( |4.1[ ) is a function with nonzero 
derivative. Since a is a polynomial in m and it is negative for m G (1/2,1) it follows that the denominator 
in (4.1) can not be zero and hence the ratio is indeed in m. A straightforward calculation shows 
that it is strictly decreasing for m G [1/2, 3/4] and strictly increasing for m G [3/4,1]. Therefore the 


proof of Theorem |2.2| applies. To complete the proof we are left to verify that the parameter intervals 
[ 5 +S’ ^ 2 fc^+ 2 /c^ ] cover the set [2/3,1). We compare the right and left endpoints of consecutive intervals 
to check if they overlap. 


fc + 2 
k + 3 ~ 


2k^ + 2fc - 1 
2 fc 2 + 2k 


0 < 2fc^ + A: - 3 


This relation holds for /c > 1, therefore all the above parameter intervals overlap. Since both endpoints 
of the fc’th interval converge to 1 as /c goes to infinity, the statement of the theorem follows. □ 
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Appendix 

Proof of Lemma \3.6\ It is enough to prove the part of the statement on the unstable direction, the stable 
part has the similar proof with T replaced by T~^. Let us first investigate the action of the tangent map 
DT on unstable vectors. At every point the unstable cone is given by the union of the first and third 
quadrants of the plane (as shown in m)- We first prove that the angle between two unstable vectors, 
lying in the same quadrant, is contracted by applying DT on the vectors. By ( |1.4[ ) it is clear that at every 
point DT has the form ( -{fab ) some quantities a, 6 < 0. The image of a unit vector (cos a, sin a) 
under the action of DT is (asina — cosa, 6 cosQ; — (1 + ab)sma). If a G ( 0 , 7 r/ 2 ) U ( 7 r, 37 r/ 2 ), then the 
original vector is an unstable vector and so by uniform hyperbolicity there exists a constant A > 1 such 
that 


(a sin a — cosa)^ + (6 cos a — (1 + ab) sina)^ > A > 1. 

The tangent map transforms the angles according to the function 

fbcosa — (1 + ab) sino' 

a —)■ arctan -^- 

V a sm a — cos a 

Differentiating this with respect to a many terms cancel out leaving 

1 


(A.l) 


(a sin a — cos a)^ + (6 cos a — (I + a 6 ) sin a)^ 


for the derivative, which is at most I/A < I for a G ( 0 , 7 r/ 2 ) U (tt, 37 r/ 2 ) using (A.l). Therefore, by the 
mean value theorem we get 
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A{D^Tvi,D^Tv2) < 


(A.2) 


for all X G A^i and vi,V 2 £ with (r>i,'y 2 ) > 0. Since we can write the difference of two unit vectors 
vi and V 2 as |tii — W 2 P = 2(1 — cosZ(z)i,X 2 )), using (A.2| we get the formula 


D^Tvi DxTv2 


\D^Tvi\ \D^Tv2\ 


1 - cos Z{DxTvi, DxTv 2 ) 
1 - cosZ(xi,X2) 


- W2p < 


< 


1 - cos(Z(ni,U 2 )/A) 
1 - cosZ(ni,t>2) 


(A.3) 


\vi - V 2 \'^ < - '*^ 2 ! 


with some constant 7 G (0,1), because ^ iTcisa^^ strictly less then 1 for every a G ( 0 , 7 r/ 2 ] and also 
in the limit as a —>■ O’*". 


We now turn (A.3) into an iterative formula for the unstable direction. Given u{x) and u{y) we estimate 
the difference of the unstable directions in T(x) and T(y), assuming that x and y are in the same 
smoothness component of the dynamics. By the invariance of unstable manifolds we know that the unit 
vector tangent to the unstable manifold at T{x) is nothing but DxTu{x) normalized. Hence we do the 
estimate 


DxTu{x) 

DyTuiy) 

< 

DxTu{x) 

DyTu{x) 


DyTu{x) 

DyTu{y) 

\DxTu{x)\ 

\DyTu{y)\ 


\DxTu{x)\ 

\DyTu{x)\ 


\DyTu{x)\ 

\DyTu{jj)\ 


(A.4) 

We remark that for this to be precise we identified the tangent spaces at x and y. Also note that the 
unstable cone field for the system is constant, so after the identification both u{x) and u{y) can be 


viewed as unstable vectors in the same tangent space. Therefore, for the second term in (A.4) we can 


apply (A.3) with x replaced by y and vi = u(x), V 2 = u{y). 


For the first term we use our assumption of x, y being in the same smoothness component and the 
assumption of Lemma |3.6| that this component is either i?o or i?i. Within these circumstances we can 
use the piecewise property of the dynamics, in particular that there exists a constant Ci > 0 such 
that 

\DxTv - DyTv\ < Cid{x, y), 


for every unit vector v. Hence for the first term of (A.4) we get 


DxTu{x) DyTu{x) 

\D^Tu{x)\ \DyTu{x)\ ~ 

^ \DxTu{x) - DyTu{x)\ \\DyTu[x)\ - \DxTu{x)\\ 
\Da:Tu{x)\ \Dx;Tu{x)\ 

< < |cid(x,y). 

|Zla;T'u(xj| A 

This, together with our previous observation on the second term leads to 


DxTu{x) 

\DxTu{x)\ 


DyTujy) 

\DyTu{y)\ 


< -Cid{x,y) + j\u{x) 


u{y)\, 


which can be iterated as long as the image of x and y are in the same smoothness component, i.e. as 
long as they are not separated by a singularity. This iteration results in 


DxTMx) 

\DxT''^u{x)\ 


DyTMy) 

iDyT-^uiy)] 




u{y)\, 


(A.5) 
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for any 0 < n < s+{x,y). To get an estimate on |M(a;) — u{y)\ (as required in Lemma (3.6)) we replace 
in (A.5) X and y by T~'^{x) and T~^{y) respectively and get 


n—1 


2=0 


|u(x) - u{y)\ < -j^Yd{T-{x),T-^{y))+rHT-{x)) - u{T-^{y))\. 

Now observe that 


(A.6) 


s+(r* '^{y)) = s+{x,y) + n- i 

S-{T^~'^{x),T^~'^{y)) = s_{x,y)+i-n, 


as long as 0 < n < S-{x, y). Hence for alH S {0,..., n — 1} we have 

s(T*"”(a;),T*"”(y)) > min{s+(x,y) + l,s_(x,?/) - n}, 
so choosing n = S-{x,y)/2 the inequality 

s(r-"(x),r-"(y))>^s(x,y) 


holds. Therefore by (3.5) still for every ie{0,...,n—l}we have 


d(r-”(x),T*-"(y)) < C9 — 


(a:,y) 


and since n > s{x,y)/2 by our choice, we get the final estimate for (A. 6 ) 

, , , , ,, 2C-\ 1 . s(3:,y) s(a:,y) 

\u{x)-u{y)\< —y^C9 2 +27 2 , 

proving the lemma with 7 ^ = max{-\/d, y^}- 


□ 
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